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Summary 

Experimental and analytical investigations focusing on aluminum/hydrocarbon gel droplet 
secondary atomization and its effects on gel-fueled rocket engine performance are being 
conducted. A single laser sheet sizing/velocimetry diagnostic technique, which should 
eliminate sizing bias in the data collection process, has been designed and constructed to 
overcome limitations of the two-color forward- scatter technique used in previous work. 
Calibration of this system is in progress and the data acquisition/validation code is being 
written. Narrow-band measurements of radiant emission, discussed in previous reports, 
will be used to determine if aluminum ignition has occurred in a gel droplet A one- 
dimensional model of a gel-fueled rocket combustion chamber, described in earlier reports, 
has been exercised in conjunction with a two-dimensional, two-phase nozzle code to 
predict the performance of an aluminum/hydrocarbon fueled engine. Estimated secondary 
atomization effects on propellant burnout distance, condensed particle radiation losses to 
the chamber walls, and nozzle two-phase flow losses also investigated. Calculations 
indicate that only modest secondary atomization is required to significantly reduce 
propellant burnout distances, aluminum oxide residual size and radiation heat losses. 
Radiation losses equal to approximately 2-13 % of the energy released during combustion 
were estimated, depending on secondary atomization intensity. A two-dimensional, two- 
phase nozzle code was employed to estimate radiation and nozzle two-phase flow effects 
on overall engine performance. Radiation losses yielded a 1 % decrease in engine I sp . 
Results also indicate that secondary atomization may have less effect on two-phase losses 
than it does on propellant burnout distance and no effect if oxide particle coagulation and 
shear induced droplet breakup govern oxide particle size. Engine 1^ was found to 
decrease from 337.4 to 293.7 seconds as gel aluminum mass loading was varied from 0- 
70 wt%. Engine I gp efficiencies, accounting for radiation and two-phase flow effects, on 
the order of 0.946 were calculated for a 60 wt% gel, assuming a fragmentation ratio of 5. 
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Program Objectives 

The overall objective of our research is to provide an increased understanding of the 
ignition and combustion characteristics of aluminum gel propellants using a combination 
of experiment and analysis, and to estimate gel combustion and secondary atomization 
effects on rocket engine performance through numerical modeling. Specific objectives are 
as follows: 

1 . Develop diagnostic techniques to investigate the ignition and secondary 
atomization of various gels in the droplet size range of practical interest (10-100 
pm) 

2. Understand the roles of surfactants, gellants, aluminum particle size, and aluminum 
mass loading on the secondary atomization of aluminum gel droplets 

3. Estimate rocket engine performance losses due to gel combustion and determine 
the effects of secondary atomization on these losses through a one-dimensional 
model of a gel-fueled rocket combustion chamber and an existing two-dimensional, 
two-phase flow nozzle code. 

4. Perform vehicle mission studies, incorporating the above losses, to estimate the 
actual performance potential of aluminum/hydrocarbon gel propellants. 

Experimental Systems 

Limitations of the Two-Color Forward-Scatter Diagnostic Technique 

In previous experimental work, a two-color, forward-scatter sizing/velocimetry system 
was used to determine particle size and velocity in a dilute stream of burning gel 
droplets. 13 He-Ne laser light scattered by a particle is detected by optics oriented at 90° 
to the laser beam axis and is used to trigger the data acquisition system, ensuring that all 
measured particles are exposed to the same incident Ar-ion light intensity, thereby 
permitting a direct correlation of particle diameter with scattered Ar-ion light intensity. 
Given the Gaussian intensity profile of the He-Ne laser beam, however, a range of trigger 
voltage levels is expected for a given size particle. Because of inherent electronic noise on 
the He-Ne trigger channel, a non-zero trigger threshold is required to prevent false 
acquisition triggering. Given that a small particle scatters less light than a large particle 
for a given incident light intensity, the effective sampling cross-sectional area in which a 
small particle can generate a sufficient trigger signal will be less than the sampling area for 
larger particles. Because the probability of a given particle size being detected is 
proportional to the sampling cross-sectional area, the small particle population is 
underrepresented. 14 Therefore, the measured particle size distribution is biased towards 
larger particles and a probability distribution function (PDF) must be applied to the data in 
order to obtain the true particle size distribution. 1 

By convolving a Gaussian beam profile and a detector slit function profile, Holve 4 
derived expressions for both sampling volume and sampling cross-sectional area as 
functions of system geometry and the light scattering response function, F(d,n), where d is 
particle diameter and n is the particle complex refractive index. Essentially, F(d,n) is the 
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ratio of scattered light intensity to incident light intensity. The detector slit function 
profile was approximated as the exponential expression, exp[-2(2x'/W s )P], where p is 
determined by the blur spot size of the receiving optics, and coordinate direction x' is at an 
angle, 9, to the coordinate x. This arrangement is depicted in Fig. 1 . 

Based on the above assumption of an exponential slit function profile, the sampling 
cross-sectional area for a particular particle size, j, can be expressed as follows: 
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where (3 is the system gain, W 0 is the He-Ne beam waist diameter, W s is the detector slit 
width, and d t is the threshold particle diameter that can be detected. This expression can 
be normalized by dividing by a reference sampling cross-sectional area, S re f, to yield the 
following: 


. P +■ 2 
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S ref ln[F(d tef ,n)/F(d t .n)] 

As indicated, Eq. 2 is the inverse of the PDF used to correct for data collection bias. 

Although Eq. 2 describes the required PDF, there are still difficulties in determining this 
expression, the greatest being uncertainty in the particle refractive index, n, and the 
scattering response function, F(d,n). Light scattering in the 90° direction, unlike near- 
forward scattering, is highly sensitive to particle refractive index and shape. Since the 
refractive index of a gel droplet is unknown and changes as the droplet passes through 
different combustion stages, the scattering response function cannot be accurately 
determined. In addition, given that the He-Ne beam diameter in the probe volume is 81 
pm, the scattering response function cannot be calculated using Mie theory because large 
droplets are exposed to a non-uniform incident illumination. Holve 4 recommends a beam 
diameter four times greater than particle diameter for uniform particle illumination, 
limiting Mie-scattering calculations for the two-color diagnostic system to particles less 
than 20 pm in diameter. 

Single Laser-Sheet Forward-Scatter Sizing/Velocimetrv Technique 

Based on the above difficulties, a single laser diagnostic technique is being developed, 
in which all particles pass through a uniform intensity region of a laser sheet, thereby 
removing any particle size bias in the data collection process. Testing has demonstrated 
that the gel droplet spray, exiting the burner through a 1 .8 mm i d. tube, maintains an 
approximately constant diameter of ~2 mm throughout the postflame region. 
Consequently, the laser sheet must have a uniform intensity region, defined as an intensity 
variation of less than 2%, large enough to encompass the entire gel particle spray. 
Calculations indicate that a 1/e 2 sheet width of ~36 mm provides a 4 mm diameter uniform 
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Figure 1 Schematic of coordinates used in Holve's calculation of sample volume 
cross-sectional area 



intensity region, which should encompass the 2 mm diameter droplet spray mentioned 
above, including particles ejected by secondary atomization. 

A schematic of the single laser sheet system and resulting laser sheet are presented in 
Figs. 2 and 3, respectively. In brief, all mm beam of a He-Ne laser (Spectra-Physics 
124B) is passed through a 750 mm focal length plano-convex spherical lens (LI, Oriel 
40815) and a 19 mm focal length plano-convex cylindrical lens (L2, Newport CKX019), 
producing a horizontal laser sheet over the burner. This sheet has a calculated 1/e 2 
thickness of 550 pm and a sheet width of 36 mm at the focal point of the 750 mm iens. 
These dimensions are depicted in Fig. 3. Gel particles passing through this sheet scatter 
light, which is collected in the near- forward direction and collimated by a 350 mm plano- 
convex spherical lens (L3, Oriel 40800) with a strip of flocking material across its center 
to block the direct laser light. This collimated light is then focused on a 200 pm horizontal 
slit (SI) by a second 350 mm lens (L4, Oriel 40800). Light passing through this slit is 
then separated into two components by a beam splitter cube (BS1, Melles-Griot 03-BSC- 
009). The first component, used for particle sizing, passes through a 632 8 nra line filter 
(FI, Oriel 52720) and enters a PMT (Hamamatsu R928) where the particle sizing signal is 
generated. The second component, used to detect aluminum combustion, passes through 
a 400 nm narrow bandpass filter (F2, Oriel 53800) and enters a second PMT (Hamamatsu 
R928) to produce the aluminum combustion signal. A two-channel 20 MHz A/D 
acquisition board (Rapid Systems 2040) will be used to sample the above signals and a 
486 33 MHz personal computer will be employed for signal validation and data analysis 
computations. 

Maximum Sample Number Density : Using the above diagnostic system, only data 
obtained from single particle light scattering provides a correct measure of particle size 
and velocity. If two or more particles are simultaneously present in the sample volume, 
the data is invalid and must be rejected. Because high rejection rates slow the acquisition 
process, a probability analysis must be performed to evaluate likely data rejection rates 
given particle number density in the flow and the sample volume size. 

Since the probability of a particle being in the sample volume at a given time is 
statistically independent of other particles being present, the probability, P(n,x), of n 
additional particles being simultaneously present when a given particle is in the probe 
volume is given by the Poisson distribution 5 
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In the above summation, x represents particle residence time in the control volume, m is 
the number of droplet size classes and C, is the count rate for the i 01 size class. Given that 
data are only valid if no additional particles are present (i.e., n=0), and that the sum of all 
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Figure 2 Schematic of the laser sheet forward-scatter diagnostic technique for droplet 
size and velocity measurements and aluminum combustion detection 







probabilities must equal unity, the probability of more than one particle being 
simultaneously present in the sample volume can be expressed as 5 


P sra , = I -exp 


- 2 > 
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The particle count rate, C p can be related to the particle number density through the 
relation, C, = U S, N,, where Li is particle velocity, S, is the probe cross-sectional area for 
particle size 1 , and N; is the i lh size class number density. Since our diagnostic technique 
employs a sample volume of uniform light intensity, S, is the same for all particle sizes 
Substituting the above relation for C p and assuming that probe volume residence time is 
independent of particle size, Eq. 5 becomes 


1-exp 
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Furthermore, the probe volume residence time, x, represents the time a particle takes to 
cross the laser sheet, t/U, where t is the sheet thickness. Using this relation and realizing 
that the probe sample volume V p =t-S, Eq. 6 can be rewritten as 

P OT1 = l-exp(-V p N T ) ( 7 ) 

where N T is the total particle number density. Figure 4 illustrates the effect of particle 
number density on the probability of n additional particles being in the sample volume. 

For n=0, no interfering particles are present in the sample volume. Therefore, the 
probability of n=0 represents the probability of obtaining a valid data point. Based on 
experience with the two-color forward scatter system, 2 3 data rejection rates should be 
acceptable for particle number densities of up to 1000 particles/cc. 1 Although this number 
density should prove adequate for our needs, experimental testing will be needed to 
provide verification. 


Theoretical Engine Performance Analysis 

One-Dimensional Combustor Model Description 

As developed in previous reports, a one-dimensional model of a JP-10/Al/O 2 (gaseous) 
fueled rocket combustor has been developed to evaluate secondary atomization effects on 
propellant combustion. 6 7 In brief, a radially uniform spray, consisting of four droplet size 
classes, enters the combustion chamber and bums in a process incorporating liquid carrier 
burnout, droplet secondary atomization, aluminum agglomerate heat up and combustion, 
two-phase particle flow, and radiation from solid combustion products to the chamber 
walls. A schematic of the gel combustion process for a single droplet, with and without 
secondary atomization, is presented in Fig. 5. Because the post-secondary atomization 
droplet size distributions are currently unknown, droplets are assumed to fragment into a 
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Figure 4. Effect of particle number density on the probability of having n 
additional particles in the probe sample volume. Results assume 
a 2 mm diameter droplet stream. 
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SECONDARY 
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and without droplet secondary atomization 


given number of equal-sized droplets where the number of secondary droplets produced 
per initial droplet is defined as the fragmentation ratio, P, which is treated as a model 
parameter. 

Radiation Modeling: A detailed development of the condensed product radiation 
model used in the combustor code is presented in prior reports, 6 - 7 but further clarification 
of this radiation model is required. In our earlier report, 7 the final equation representing 
radiation heat flux from the condensed products to the chamber walls was expressed as 


q r (R) = Z 




( 8 ) 


where n represents a wavelength interval; I 0 and I,, Modified Bessel functions; £,, a 
function of particle optical properties; T the gas temperature; T w , the chamber wall 
temperature; and F n (T), the fraction of total radiation emitted by a blackbody at 
temperature T in the wavelength interval between X n and X n+1 . The variable R was defined 
as the chamber radius. It should be clarified that R is actually the chamber radius in terms 
of flow optical thickness and should be defined as R =(a + a s )R, where R is the physical 

chamber radius and the variables, a and ct s , are the absorption and scattering coefficients of 
the flow. 

One-Dimensional Combustor Model Results 

The one-dimensional combustor code was exercised using the conditions presented in 
Table 1. These values were chosen to simulate an upper-stage booster and represent the 
maximum I sp operating point for a 60 wt% aluminum gel. Because propellant mass flow 
rate is governed in part by nozzle geometry, the Solid Propellant Performance Prediction 
Code (SPP) 8 was used to determine the propellant mass flow rates. Since an appropriate 
spray distribution remains to be determined from gel atomization research, the arbitrary 
normalized droplet size distribution presented in Fig. 6 was used. 


| Table 1: Model Operating Conditions 

Chamber Diameter 

0.19 m 

Chamber Pressure 

38 atm 

Liquid Carrier 

JP-10 

Gel Aluminum Mass Loading 

60% 

Gel Flow Rate 

10.78 kg/s 

Oxidizer Flow Rate 

1 1.75 kg/s 

Chamber Wall Temperature 

1000 K 


Gas temperature and composition, assuming a fragmentation ratio, P, of 5 are 
presented in Fig. 7. The jagged shape of both the temperature and composition profiles 
arise from using only four droplet size classes instead of a continuous droplet size 
distribution. The slow rate of temperature increase in the region between 0. 1 and 0 25 m 
is caused by the large enthalpy transfer from the gas flow to heat the aluminum after 
hydrocarbon burnout. 
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NORMALIZED NUMBER OF DROPLETS 



DROPLET SIZE (jim) 


Figure 6. Normalized droplet size distribution used in the combustor code. 
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AXIAL LOCATION (m) 


Figure 7. Gas temperature and major species mass fluxes versus axial 
location. Data are for a 60 wt% aluminum gel assuming a 
fragmentation ratio of 5 






In Fig. 8, propellant burnout distance is plotted versus fragmentation ratio, [3, to illustrate 
the potential benefits of secondary atomization in reducing required engine residence 
times. It is readily apparent that only modest secondary atomization is required to 
significantly reduce propellant burnout distance, and that higher secondary atomization 
intensities, represented by larger fragmentation ratios, have a marginal effect 

Final oxide residual diameter as a function of fragmentation ratio is shown in Fig 9 
Similar to the trend seen in propellant burnout distance, moderate secondary atomization 
intensities significantly reduce final A1 2 0 3 residual diameter, with greater secondary 
atomization intensities providing decreasing marginal reductions in residual diameter The 
final A1 2 0 3 weight distribution as a function of particle size is plotted in Fig. 10. The bi- 
modal size distribution associated with aluminum combustion 910 is readily apparent as well 
as the decrease in residual particle size caused by secondary atomization. It should be 
noted that secondary atomization has no effect on the size of the sub-micron fume 
particles produced through vapor-phase oxidation of the aluminum. 

Radiation losses from condensed combustion products to the chamber walls as 
functions of fragmentation ratio and aluminum mass loading are shown in Fig. 1 1 . 

Chamber operating conditions are those that produce maximum 1^ for each aluminum 
mass loading. Radiation losses were found to be strongly influenced by the fragmentation 
ratio, primarily due to changes in propellant burnout distance. Increasing the gel 
aluminum mass loading for a given fragmentation ratio was also found to increase 
radiation losses because of increased chamber temperature and slight increases in 
propellant burnout distance with increasing aluminum mass loading. Radiation losses 
were approximately equal to 2- 1 3 % of the energy released during combustion, depending 
on the aluminum mass loading and the fragmentation ratio. 

It should be cautioned that the above radiation heat transfer calculations are only 
approximate. Axial radiation heat transfer, which has been neglected, may alter the 
droplet combustion process and/or the actual radiation losses. In addition, research 11 has 
indicated that the optical properties of the condensed combustion products can vary 
considerable depending on propellant composition and oxidizer/ftiel mixture ratio 

SPP Modeling of Nozzle Two-Phase Flow and Engine Performance 

Two-phase flow effects in the engine nozzle were estimated using the SPP nozzle 
code, 8 which contains two solution modules for determining nozzle performance. The 
first module provides a transonic analysis of the converging and throat region of the 
nozzle, and the second module provides a Method of Characteristics (MOC) solution of 
the supersonic portion of the nozzle. This MOC module incorporates the following 
assumptions: 


• Gas and particle properties are variable 

• Energy transfer between the gas and particles is only through convection 

• Oxide particles are spherical 

• Particles do not interact 

• Mass is not transferred between the particles and the gas 

• Condensed phase occupies negligible volume in the flow 


13 




5 



FRAGMENTATION RATIO 


Figure 8. Propellant burnout distance as a function of secondary atomization 
fragmentation ratio for a 60 wt% aluminum gel. 
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Figure 9. Final oxide residual diameter as a function of secondary atomization 
ratio for a 60 wt% aluminum gel. 
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DIAMETER (|im) 


Figure 10 Aluminum oxide weight fraction versus particle diameter for a 
60 wt% aluminum gel. 
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Figure 1 1 Radiation heat loss as a percentage of combustion energy released 
as a function of fragmentation ratio and gel aluminum mass loading. 
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It should be noted, however, that because particle interaction and mass transfer between 
the particles and gas are neglected in the SPP code, the A1 2 0 3 residual size distribution 
predicted by the combustor code may not be correct for nozzle performance calculations 
Since the A1 2 0 3 particles are molten throughout most of the nozzle, and because small 
particles accelerate more quickly than large particles, particle size may increase through 
coagulation. 12 Similarly, additional A1 2 0 3 may be produced through the recombination of 
gas-phase radicals as the exhaust gases cool during expansion, resulting in the nucleation 
of additional particles and/or growth of previously formed particles. 13 Particle size may 
also decrease due to shear breakup of droplets, particularly in the throat region of the 
nozzle. 14 

Because of these uncertainties in A1 2 0 3 particle size, two methods of estimating 
particle size, which should bound the true particle size, are used in the evaluation of nozzle 
performance. In the first method, we use an A1 2 0 3 particle size distribution that is 
determined by the one-dimensional combustor code, making secondary atomization the 
primary mechanism governing particle size In the second method, secondary atomization 
is assumed to have no effect on mean particle size; rather, coagulation, particle surface 
growth, and shear induced droplet breakup are assumed to be the dominate mechanisms 
affecting particle size. Assuming that these coagulation and breakup mechanisms are 
comparable to those in solid nozzles, the following solid motor correlation of mass median 
A1 2 0 3 particle diameter, D 43 , as a function of nozzle diameter, D t , was employed: 8 

D 43 = 3.63D, 2932 (9) 

resulting in a particle mass median diameter of 5.6 pm. 

Engine Performance Results 

A comparison of two-phase flow effects on 1^ for the above methods of determining 
A1 2 0 3 particle size are presented in Fig. 12 for a 60 wt% aluminum gel. As seen here, 
secondary atomization may reduce two-phase flow losses but not as significantly as 
propellant burnout distance (4% versus 62%, cf. Fig. 3). To illustrate the separate 
contributions of radiation and two-phase losses, radiation losses were incorporated in the 
case represented by the dashed line. Given a fragmentation ratio of 5, radiation losses 
yield a decrease in 1^ of approximately 1 % compared to the 4% resulting from two-phase 
losses. 

Using the solid motor correlation for particle size, 1^ was calculated for a range of 
aluminum mass loadings and propellant mixture ratios and compared with I sp calculations 
for a JP-10/O 2 bi-propellant. Figure 13 shows that 1^ decreases with increasing aluminum 
mass loading and that the maximum I sp mixture ratio becomes richer, as has been predicted 
by other studies, 1516 although the 1^ decrease is significantly greater with the 
incorporation of radiation and two-phase flow losses. 
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0/F MIXTURE RATIO 


Figure 12. Engine I sp as a function of oxide particle size model and 
O/F mixture ratio for a 60 wt% aluminum gel. Dashed 
line indicates inclusion of radiation losses. Model assumes 
an Extended Delta nozzle profile. 
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0/F MIXTURE RATIO 


Figure 13. Engine I sp as a function of O/F mixture ratio and gel 

aluminum mass loading. Dashed lines represent inclusion 
of radiation losses. 




Conclusions 


Based on the above experimental and theoretical work, the following results were 
obtained: 

1 A single laser-sheet sizing/velocimetry technique has been developed to overcome 
limitations in the two-color, forward-scatter diagnostic technique used in previous 
work. This system has been constructed and is now being calibrated. Data 
acquisition/validation software is currently being written to operate this system. 

2. Exercise of the one-dimensional combustor model predicts that only moderate 
secondary atomization ((3=5) is required to reduce overall propellant burnout 
distance by 62% and final A1 2 0 3 residual diameter by 41%. 

3. Radiation losses for a 60 wt% gel, assuming (3=5, are estimated to be approximately 
5 % of the energy released during combustion, resulting in a 1% decrease in engine 
I sp . Secondary atomization may reduce radiation heat transfer losses, primarily due 
to decreases in the propellant burnout distance. For example, a fragmentation ratio 
of five results in a 6 1 % decrease in radiation losses. Radiation losses were also 
found to be a function of gel composition and engine operating conditions, as a 
result of changes in chamber temperature. 

4. Two-dimensional, two-phase nozzle code results indicate that secondary 
atomization may have little effect on nozzle two-phase flow losses. Specifically, a 
fragmentation ratio of 5 decreases two-phase flow losses by only 4 % compared to 
the above 62% reduction in propellant burnout distance. Furthermore, secondary 
atomization may have no effect on two-phase flow losses if particle coagulation, 
surface growth, and shear induced breakup are the dominant mechanisms controlling 
oxide particle size. 
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Future Work 


1 . Calibrate single-laser diagnostic apparatus and complete the data acquisition code. 

2. Perform a parametric study of aluminum mass loading and surfactant concentration 
effects on limiting droplet size for secondary atomization and post-secondary 
atomization particle size distributions. 

3. Numerically model the burning droplet spray used in the above experiments using a 
modified version of the one-dimensional combustor code to aid in differentiating 
between normal droplet combustion and secondary atomization. 

4. Using data collected in the above work, incorporate a realistic post-secondary 
atomization model in the one-dimensional combustor code and using the 
combustor code and SPP, evaluate the potential of gel-fueled engine performance 
and secondary atomization effects on this performance. 

5. Develop a realistic Earth-to-LEO mission analysis to estimate the performance 
potential of hydrocarbon gel fueled vehicles. 
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